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Steep gradients of temperature and density, called cold fronts, are observed by 
Chandra in a leading edge of subclusters moving through the intracluster medium 
(ICM). The presence of cold fronts indicates that thermal conduction across the 
front is suppressed by magnetic fields. We carried out three-dimensional mag- 
netohydrodynamic (MHD) simulations including anisotropic thermal conduction 
of a subcluster moving through a magnetically turbulent ICM. We found that 
turbulent magnetic fields are stretched and amplified by shear flows along the 
interface between the subcluster and the ambient ICM. Since magnetic fields re- 
duce the efficiency of thermal conduction across the front, the cold front survives 
at least 1 Gyr. We also found that a moving subcluster works as an amplifier 
of magnetic fields. Numerical results indicate that stretched turbulent magnetic 
fields accumulate behind the subcluster and are further amplified by vortex mo- 
tions. The moving subcluster creates a long tail of ordered magnetic fields, in 
which the magnetic field strength attains (3 = -P gas /-Pmag ^ 10. 
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Subject headings: conduction — MHD — galaxies: magnetic fields — X-rays: galax- 
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Introduction 



X-ray observations of clusters of galaxies by ASCA revealed complex temperature dis- 
tributions in the intracfuster medium (ICM), in wh ich hot and cool plasma coexists (e.g, 
Perseus Cluster, lArnaud et al.l 1 1994 iFurushol l200ll ). Sharp discontinuities of density and 
temperature in ICM, called cold fronts, were found by high spatial resolution observations 
by Chandra ( IMarkevitch et al.ll2000l : IVikhlinin et al.l l2001bl ). Cold fronts manifest the co- 
existence of hot and cool plasma in clusters of galaxies. They provide a key to understand 
thermal properties of the ICM. 

Cold fronts in merging clusters such as A2142, A3667, and 1E0657-56 result from merg- 
ing. When a subcluster is moving in the ICM, a sharp boundary is formed between the 
subcluster and the ambient hot ICM in the forehead of the subcluster because the cold 
plasma confined by the subcluster is subjected to the ram pressure. The cold fronts are not 
shock fronts because the Chandra images of the X-ray surface brightness show that the tem- 
perature decreases on the denser side. A3667 has a clear, large-scale (~ 500 kpc), arc-shaped 
cold front which shows a steep gradient of the X-ray surface brightness and temperature. 
The temperature decreases toward the denser part from 8keV to 4keV wit hin 5 kpc. This 



thickn ess of the front is 2-3 times smaller than the Coulomb mean free path (IVikhlinin et al. 
2001bh . 



A question which needs to be answered is how the temperature gradient is created and 
sustained in the ICM which typically has high Spit zer conductivity, 

1 K -1 (|Spitzerlll962l). Thermal conduction rapidly smooths 



«Sp = 5 x 10 7 T 5//2 ergs 1 cm 



such a steep gradient (e.g., iTakahara &: Ikeuchil Il977l ). The time required for heat to dif- 
fuse by conduction across a length L in the ICM is roughly given by Ts p ~ pL 2 / ks p 



10 7 (fcr /5k eV)~ 5/2 (n/10~ 3 cm~ 3 )(L/ 100 kpc) 2 yr, where p is the density. Ettori k Fabian 



(120001 ) and IMarkevitch et al.l (120031 ) estimated that the effective thermal conduction is at 
least an order of magnitude lower than the Spitzer value. This suggests that t he thermal con 



duction across the front i s suppr essed by magnetic fields parallel to the front (IVikhlinin et al. 



2001bl ). IVikhlinin et al.l (j2001al ) suggested that ordered magnetic fields are formed in front 
of the subcluster because small-scale turbulent magnetic fields are compressed and stretched 
along the front ahead of the su bcluster by its motion. Typical clusters o f galaxies possess 



magne tic fields of ~ fiG (e.g., iKronbergl Il994l : ICarilli fc Taylorl 120021 ) . iJohnston-Hollitt 



(120041 ) reported that 1 — 2 fiG fields, tangled on 100 kpc scale pervade the central region of 
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A3667. When magnetic fields exist, the characteristic scale of the heat exchange across the 
field lines is reduced significantly compared with non-magnetized ICM to the Larmor radius, 
?"l ~ 2500 /iG) _1 (T/5 keV) 1,/2 km. Intracluster magnetic fields play a crucial role for 
thermal conduction even if the magnetic pressure is lower than the gas pressure. 

A number of authors reported t hat cold fronts were reproduced in numerica l simulations 
as a result of a merging process (e.g.. iBialek et al.ll2002t iNagai fc Kravtsovll2003l ; iHeinz et al. 
20031 ; lAcreman et al.ll2003l ; iTakizawal 120051 ) . In these simulations, however, magnetic fields 
and thermal conduction we r e not included. To study the evolution of intracluster mag- 
netic fields, iRoettiger et al.l (119991 ) performed MHD simulations of merging clusters and 
Dolag et al. ( 20021 ) p erformed c o smolo gical MHD simulations. However, they ignored the 
thermal conduction. lAsai et al.l (120041 ) performed two-dimensional (2D) MHD simulations 
of a subcluster moving through uniform magnetic fields by including anisotropic thermal 
conduction. They showed that ordered magnetic fields wrap the subcluster and suppress 
the thermal cond uction across the front. This work was extended to three-dimension by 
Asai et all J2005h . 



Some authors pointed out that the effective conductivity in turbulent ma g netic fields is 
only s everal times lower than the Spitzer value (INarayan fc Medvedevl l200ll ). iDolag et al. 
(120041 ) carried out cosmological hydrodynamic simulations including the thermal conduction 
with the isotropic effective conductivity, k ~ and showed that temperature gradients 

are smoothed compared to the case without thermal conduction. However, a moving subclus- 
ter may stretch the turbulent magnetic fields and such fields may reduce the conductivity. 

Lyutikov (2006) studied the magnetic draping mechanism by which a strongly mag- 
netized, thin boundary layer with a tangential magnetic field is created in the boundary 
between a moving cloud and the ambient plasma. He theoretically predicted that for super- 
sonic cloud motion, magnetic field strength inside the layer reaches near equipartition values 
with thermal pressure. 



Asai et al.l (120041 ) carried out a simulation starting from a disordered magnetic field and 



showed that cold fronts can be sustained because magnetic fields are stretched along the 
front and wrap the subcluster. However, their simulation was limited to 2D and their initial 
magnetic field still had a large coherent length. Thus magnetic fields can easily suppress the 
thermal conduction across the front. 

In this paper, we present results of three-dimensional (3D) MHD simulations of a sub- 
cluster moving through turbulent magnetic fields to show that even when magnetic fields are 
turbulent, cold fronts can be sustained. The paper is organized as follows. In §2], we describe 
the initial condition and model parameters. Simulation results are presented in $3j we show 
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the effects of magnetic fields on cold fronts and amplification of magnetic fields behind the 
subcluster. Finally, we discuss and summarize our results in §H 



Simulation Model 



We simulated the time evolution of a cluster plasma in a frame comoving with the 
subcluster. The basic equations are as follows: 



dp 
dt 



+ V -(pv) = 0, 



dv 



V) v 
dB 



-Vp + 



fV x B) x B 



dt 



V x ( v x B) 
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dt 



- 2 pv 
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(1) 
(2) 

(3) 

-pv ■ V^, 

(4) 



where p, t>, p, £?, and ip are the density, velocity, pressure, magnetic fields, and gravitational 
potential, respectively. We use the specific heat ratio 7 = 5/3. The subscript || denotes 
the components parallel to the magnetic field lines. We assume that heat is conducted only 
along the magnetic field lines. We solved equations (CO)-® in a Cartesian coordinate system 
(x, jy , z) by using a solver ba sed on a modified Lax-We ndroff method (IRubin fc Burstein 
19671 ) with artificial viscosity ( Richtmyer fc Morton 1967 ) implemented to the Coordinated 
Astronomical Numerical Software (CANS). We did not include the physical viscosity. Arti- 
ficial viscosity is included only in regions close to the discontinuities to suppress numerical 
oscillations. It does not affect the dynamics in smooth regions. The thermal conduction 
term in the e nergy equation is solved by the implicit red and black successive over-relaxation 



method (see lYokoyama fc Shibatal l200ll . for detail). The radiative cooling term is not in- 



cluded. The units of length, velocity, density, pressure, temperature, and time in our simu- 



lations are r = 250 kpc, v ■ 
kT = 4keV, and t = r /v 



1000 km s 1 , p = 5 x 10 27 gem 
3 x 10 8 yr, respectively. 



3x10 11 erg cm 



We carried out simulations for 6 models. Table [T] shows the model parameters. When 
magnetic fields exist, heat conducts only parallel to the magnetic field lines. Model MT1 and 
MT2 are models with turbulent magnetic fields. We define Fourier components of magnetic 
vector potential, A(k) = Aok~ a , where k is a wave number. The amplitudes Aq are taken 
to be random by using random numbers, and we adopt a = 5/3. This vector potential A in 
/c-space is transformed to vector potential A in physical space via a 3D fast Fourier transform 
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(FFT). We computed a tangled divergence-free initial magnetic field via B = V X A. On 
the other hand, models MU1 and MU2 are models with uniform magnetic fields parallel to 
the ^-direction and perpendicular to the motion of the subcluster. 

Figure [1] shows the initial density distribution for model MT1 at z = plane. Solid 
curves and arrows show the contours of magnetic field strength and velocity vectors, re- 
spectively. We assume that a spherical isothermal low-temperature (kT in = 4keV) plasma 
is confined by the gravitational potential of the subcluster. The subcluster has a /3-model 
density distribution, 



where we adopted f3' = 2/3, the core radius r c = 290 kpc, and the maximum density p c = 
2 p = 10~ 26 gcm~ 3 . The subcluster is embedded in the low density (p ou t = Po/4), hot 
(kT out = 2 kT in ) ambient plasma. We assume that the subcluster is initially in hydrostatic 
equilibrium and has a jump of density and temperature with respect to the ambient plasma. 
We also assume that the ambient plasma has a uniform speed M = v x /c St0nt = 1, where 
c SjOUt is the ambient sound speed. Note that magnetic fields exist even inside the subcluster 
in all models with magnetic fields. The box size of our simulations is (2.5 Mpc) 3 . We used 
256 3 grid points for typical models. The numerical resolution is ~ 10 kpc. 

An important parameter is plasma (3 defined as the ratio of the gas pressure to the 
magnetic pressure (see Tabled]). The initial mean field strength Bo is ~ 0.03 pG for model 
MT1 and ~ 0.09 pG for model MT2. Correspondingly, initial plasma (3 for these models is 
(3q ~ 7.5 x 10 4 and ~ 8.3 x 10 3 , respectively. For models with uniform fields, the strength 
B is ~ 0.07 pG for models MU1 and ~ 0.27 pG for model MU2. The initial plasma /3 for 
these models is (3q ~ 1.0 x 10 4 and ~ 1.0 x 10 3 , respectively. Model MT3 is the same as 
model MT1 except that thermal conduction is ignored and the box size of the x-direction 
is 1.5 times larger than that of model MT1. We also carried out a simulation for a model 
without magnetic fields (model H), including isotropic thermal conduction. 

For turbulent field models (MT1, MT2, and MT3), the left boundary at x — —5 is 
taken to be the fixed boundary, except for the magnetic fields. The magnetic fields at the 
left boundary are extracted from the initial distribution of magnetic fields as follows, 




(5) 



1, j, k) = B (ix - v x0 t, j, k), 



(6) 



where B is the initial magnetic field, j and k are mesh numbers in the y and z directions, 
ix is the total number of grid points in the x-direction, v x q is the initial velocity of the 
subcluster, and t is the time. For other models, the left boundary is taken to be a fixed 
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boundary. In all models, boundaries other than the left boundary are free boundaries where 
waves can be transmitted. 



3. Results 

3.1. Time Evolution of a Subcluster and Magnetic Fields 

Figure [2] shows time evolution of density distribution and magnetic fields for model 
MT1. The left panel shows the initial state and the right panel shows the distribution at 
t = 1.0 Gyr. Solid curves show the magnetic field lines. Since the subcluster plasma moves 
with sound speed, a bow shock appears ahead of the subcluster. Magnetic field lines are 
stretched along the subcluster surface at t — 1.0 Gyr due to the ambient gas motion. Thus, 
the motion of the subcluster creates ordered magnetic fields along the interface between the 
subcluster and ambient ICM. The magnetic fields accumulate behind the subcluster. 

3.2. Effect of Magnetic Fields on Thermal Conduction 

Figure [3] shows snap shots of distributions of temperature (left panels), temperature 
gradients (middle panels), and magnetic field strength (right panels) at t — 1.0 Gyr. The 
upper panels show the slices at z = plane and the bottom panels show the slices at 
x = plane for model MT1, respectively. Solid curves in the left panels are the contours of 
magnetic field strength. Arrows in the left (and right) panels show the velocity vectors, and 
those in the middle panels are the gradients of temperature. 

The temperature distributions (left panels) show that steep temperature gradients are 
maintained at t — 1.0 Gyr because the thermal conduction across the front is suppressed by 
stretched magnetic field lines wrapping the subcluster. The middle panels show that steep 
temperature gradients around the subcluster surface are sustained. The cold front is located 
at (x, y) ~ (— 1, 0) in z = plane (upper panels). We can also identify the bow shock at 
(x, y) ~ (—3, 0) in z = plane. The steepest temperature gradient is seen in the vicinity 
of (x, y) = (—1, 0) in the upper middle panel. The distributions in x = plane (bottom 
panels) clearly show that the subcluster is almost entirely covered with magnetic fields. 

In the left and right panels of Figure [3, we can see that magnetic fields are stretched 
and compressed around the subcluster as we already mentioned. A shear flow along the 
boundary between the subcluster and the ambient plasma stretches magnetic fields. More- 
over, magnetic field strength is amplified behind the subcluster (see §3.31 for details). The 
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field amplification is more prominent in the tail of the subcluster than in the forehead. The 
amplification of magnetic fields ahead of the cold fronts is due to the compression of magnetic 
fields by the ambient plasma flow hitting the subcluster. The ambient plasma flowing along 
the subcluster surface converges to the x-axis behind the subcluster. Therefore, magnetic 
fields frozen to the plasma accumulates behind the subcluster. 

In order to demonstrate the effects of magnetic fields on thermal conduction, we present 
results of hydrodynamic model (model H) in Figure HI The left and right panels show the 
distributions of temperature and temperature gradients. Compared with the results shown in 
Figure [3X temperature gradients are smeared out because isotropic thermal conduction from 
the ambient hot plasma rapidly heats up the dense cool plasma confined in the subcluster. 

In Figure 0, we show the distributions of quantities at t = 1.0 Gyr along the x-axis 
(y — z — 0). The left and right panels show the distributions of temperature (solid curve), 
density (dashed curve), and pressure (dotted curve) for models MT1 and H. The temperature 
distribution in the left panel shows that a steep gradient exists at x = — 1, while pressure 
distri bution is smooth. This feature is consistent with the observed features of cold fronts 
(e.g., iMarkevitch et al.ll2000l ). On the other hand, when magnetic fields do not exist (right 
panel), the subcluster plasma is subjected to the isotropic thermal conduction. After t = 
1.0 Gyr, the subcluster evaporates because of the conduction from the ambient hot plasma. 
The peak density in the subcluster becomes lower than that in the initial state. 



3.3. Amplification of Magnetic Fields 

Simulation results revealed other interesting features. We found that a moving sub- 
cluster works as an amplifier of magnetic fields. The left panel of Figure [6] shows the 
distributions of magnetic field strength for the turbulent field models (models MT1 and 
MT2). The right panel shows the distributions of plasma (3. We plot these distributions 
along (y, z) = (-0.04, -0.20) for model MT1 (black curves) and (y, z) = (-0.04, -0.24) 
for model MT2 (blue curves), respectively. Solid curves in both panels are distributions at 
t = 1.0 Gyr and dashed curves show those at the initial state, respectively. When turbulent 
fields exist, field strength is amplified in front of the subcluster and behind it. The amplifi- 
cation of the field strength is most prominent behind the subcluster. The field strength in 
both models increases about 30 times with respect to the averaged initial value. The right 
panel shows that plasma j3 in both models decreases. In model MT2, plasma j3 decreases 
below P ~ 10 in the tail of the subcluster. 

Let us compare the results for the turbulent field models with those of uniform field 
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models. Figure [7] shows the distribution of plasma j3 at t — 1.0 Gyr along y = plane 
for model MT2 (left) and that for model MU2 (right). Note that the initial direction of 
magnetic fields for model MU2 is parallel to the z-axis. In both panels, plasma (3 decreases 
remarkably behind the subcluster. The plasma f3 in this region is lower than f3 ~ 10. The 
region of lower plasma (3 is larger for model MU2 than that for model MT2 because the 
initial magnetic fields are uniform in model MU2, thus the ambient plasma flow creates the 
ordered fields easily behind the subcluster. An important finding is that plasma (3 decreases 
to (3 ~ 10 behind the subcluster even if the initial magnetic fields are turbulent. 

In model MT3, we used larger simulation box and carried out a simulation until t = 
3.0 Gyr in order to follow the growth of magnetic fields behind the subcluster. The size 
of simulation box for model MT3 is 1.5 times larger in the x-direction than that for other 
models. Figure [S] shows the distributions of magnetic field strength (left) and plasma f3 
(right) along the x-direction (y = z = —0.08) for model MT3. Curves in both panels show 
the distribution at t = 0.0 Gyr (thin solid curve), 1.0 Gyr (dashed curve), 2.0 Gyr (dotted 
curve), and 3.0 Gyr (thick solid curve). Magnetic fields are amplified with time behind the 
subcluster. They are accumulated around x ~ 3 at t = 1.0 Gyr. The flow motion stretches 
magnetic fields along the x-direction. At t = 3.0 Gyr, magnetic field strength peaks around 
x ~ 7. This strength is about 80 times higher than the initial value. In the right panel, the 
minimum plasma f3 appears at x ~ 6. In this region, the plasma j3 is smaller than 10. 

Figure M shows the distributions of plasma (3, velocity vectors, and magnetic field vectors 
in z = plane at t — 3.0 Gyr for model MT3. The region behind the subcluster, 2.5 < x < 10, 
—2.5 < y < 2.5, z = is shown. The front of the subcluster is not shown here (it is located 
at x ~ 0). The top panel shows that (3 < 10 along the x-axis (y = z — 0). We found that 
the ambient gas passing through the subcluster generates vortices behind the subcluster. 
They create back flows a l ong t he x-axis (see middle panel). The back flow is similar to that 



reported by iHeinz et al.l (120031 ) . Magnetic fields dragged by the gas flow are thus stretched 
along the x-axis behind the subcluster. The x-component of magnetic fields is particularly 
strengthened in the vicinity of the x-axis (see bottom panels). 



4. Discussion and Summary 

We carried out 3D MHD simulations of a subcluster moving through a magnetically 
turbulent ICM in order to study the effects of magnetic fields on thermal conduction. We 
assumed that a dense, cold subcluster with a sharp discontinuity in density and temperature 
is moving with sound speed. 
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In this paper, we studied whether a cold front can be sustained in a magnetically 
turbulent ICM. In §3.11 and §3.21 we demonstrated that a cold front is maintained for over 
1 Gyr because magnetic fields stretched along the front suppress the thermal conduction 
across the front even if magnetic fields are initially turbulent. On the other hand, when 
magnetic fields do not exist, steep temperature gradients cannot be maintained because the 
cold subcluster is heated up due to the isotropic thermal conduction. Therefore, we conclude 
that cold fronts in merging clusters can exist because magnetic fields coupled with the motion 
of the subcluster suppress the thermal conduction. 

In §3.3[ we showed that magnetic fields are amplified significantly behind the subcluster 
because the ambient flow converges to the tail of the subcluster. The flow accumulates 
magnetic fields to this tail region. Furthermore, the vortex motions behind the subcluster 
accumulates and stretches the magnetic fields. The enhanced fields are maintained for a 
long time behind the subcluster. Thus, the motion of a subcluster forms a long tail of 
ordered magnetic fields. It may be worth noting that the magnetic filaments created behind 
the subcluster are similar to the magnetic filaments in the solar atmospheres, in which the 
magnetic fields are accumulated in boundaries of the convective cell. Plasma (3 decreases 
to (3 < 10. If such a long tail of magnetic fields interacts with other moving subclumps, 
magnetic fields will be further amplified. Consequently, small-scale weak magnetic fields in 
the ICM can be amplified and create large-scale magnetic fields whose energy is comparable 
to the thermal energy. This mechanism may also apply to the amplification of small-scale 
primordial magnetic fields once dark matter clumps are formed. 

We thank T. Yokoyama for developments of the coordinated astronomical numerical 
software (CANS) which include 2D and 3D MHD codes including thermal conduction. The 
development of CANS was supported by ACT-JST of Japan Science and Technology Corpo- 
ration. This work is supported by JSPS Research Fellowships for Young Scientists. Numerical 
computations were carried out on VPP5000 at the Center for Computational Astrophysics, 
CfCA, of the National Astronomical Observatory of Japan and joint research program of 
IMIT, Chiba University. 
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Table 1: Simulation models and parameters. 



Model 


K 3 - 


magnetic field 


5 [/xG] b 


A) c 


box size [Mpc 3 


number of grids 


MT1 


K\\ 


turbulent 


0.03 


7.5 x 10 4 


2.5 3 


256 3 


MT2 


K\\ 


turbulent 


0.09 


8.3 x 10 3 


2.5 3 


256 3 


MU1 


K\\ 


uniform 


0.07 


1.0 x 10 4 


2.5 3 


256 3 


MU2 


K\\ 


uniform 


0.27 


1.0 x 10 3 


2.5 3 


256 3 


MT3 





turbulent 


0.03 


7.5 x 10 4 


3.75 x 2.5 x 2.5 384 x 256 x 256 


H 


K 







oo 


2.5 3 


256 3 



a K is the thermal conductivity, and the subscript || denotes the component parallel to magnetic field lines. 
b ' c Bg and (3q are initial mean magnetic field strength and initial mean plasma (3. 
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t=0.0 log p 




Fig. 1. — Initial distribution of logarithm of density at z = plane. Solid curves and arrows 
show the contours of strength of magnetic fields and velocity vectors, respectively. 



f = 0.0 Gyr , t= 1.0 Gyr 




Fig. 2. — Time evolution of density distribution and magnetic field lines. Left and right 
panels show the distributions of density and magnetic field lines at t = 0.0 Gyr and t = 
1.0 Gyr, respectively. Color shows isosurfaces of density (logp). Solid curves show magnetic 
field lines. 
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Fig. 3. — Distributions of temperature (left panels), gradients of temperature (middle pan- 
els), and magnetic field strength (right panels) at t = 1.0 Gyr for model MT1. Upper panels 
show the slice at z = and lower panels show the slice at x = 0. Solid curves in left panels 
show the contours of magnetic field strength. Arrows in the left and right panels show the 
velocity vectors. Arrows in the middle panels show the gradients of temperature. 
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Fig. 4. — Distributions of temperature (left) and gradients of temperature (right) in z = 
plane at t — 1.0 Gyr in hydrodynamic model (model H). Arrows show the velocity vectors 
(left) and the gradients of temperature (right), respectively. 
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Fig. 5. — Comparison of distributions of temperature, density, and pressure along x-axis 
(y = z = 0) at t = 1.0 Gyr for models MT1 (left) and H (right). Solid, dashed, and dotted 
curves show the temperature, density, and pressure, respectively. 




Fig. 6. — Distributions of magnetic field strength (left) and plasma (3 (right) for models 
MT1 and MT2. Black curves in both panels show the distributions along the x-direction, 
(y, z) = (—0.04, —0.20) for model MT1. Blue curves in both panels show the distributions 
along the x-direction, (y, z) = (—0.04, —0.24) for model MT2. Solid and dashed curves in 
both panels show those at t — 1.0 Gyr and the initial state, respectively. 
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Fig. 7. — Distributions of plasma (3 in y — plane at t — 1.0 Gyr for model MT2 (left) and 
model MU2 (right). Arrows show the velocity vectors. 




Fig. 8. — Distributions of magnetic field strength (left) and plasma f3 (right) for model MT3 
along x-direction. We plot them along y = z = —0.08 for both panels. The length of the 
x-axis is 3.75 Mpc. Curves show the distributions for t = 0.0 Gyr (thin solid curve), 1.0 Gyr 
(dashed curve), 2.0 Gyr (dotted curve), and 3.0 Gyr (thick solid curve), respectively. 
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Fig. 9. — Distributions of plasma /?, velocity vectors, and magnetic field vectors from top 
to bottom in z = plane at t = 3.0 Gyr for model MT3. The region behind the subcluster 
[x, y, z] = [2.5 : 10, —2.5 : 2.5, 0] is shown. The front of the subcluster around x ~ is not 
shown here. Arrows in middle and bottom panels show velocity vectors and magnetic field 
vectors. 



